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PREFACE 


This  report  documents  the  General  Analysis  Module  of  the  three- 
dimensional  stability  analysis/design  (3DSAD)  program.  The  module  was 
developed  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES), 
Vicksburg,  Miss.,  in  the  Automatic  Data  Processing  (ADP)  Center  by  Mr.  Fred  T. 
Tracy.  The  work  was  sponsored  through  funds  provided  WES  by  the  Civil  Works 
Directorate,  Office,  Chief  of  Engineers,  U.  S.  Army  (OCE),  under  the  Computer- 
Aided  Structural  Engineering  (CASE)  Project. 

Specifications  for  the  program  were  provided  by  the  members  of  the  CASE 
Task  Group  on  3D  Stability.  The  members  of  the  task  group  during  the  period 
of  development  were  as  follows: 

Mr.  Charles  W.  Kllng,  Mobile  District  (Chairman) 

Mr.  Robert  Haavisto,  Sacramento  District 
Mr.  John  Hoffmeister,  Nashville  District 
Mr.  Gerrett  L.  Johnson,  Seattle  District 
Mr.  Thomas  J.  Mudd,  St.  Louis  District 
Mr.  William  Holtham,  New  England  Division 

This  report  was  written  by  Mr.  Tracy  and  Mr.  Charles  W.  Kllng,  Mobile 
District.  Mr.  Kllng  also  contributed  significantly  as  an  experienced 
practicing  engineer.  Mr.  Donald  R.  Dressier,  Structures  Division,  Civil  Works 
Directorate,  OCE,  and  Dr.  Howard  B.  Wilson,  Jr.,  University  of  Alabama,  are 
acknowledged  for  their  contributions  to  the  kern  plot.  Ms.  Debra  G. 

Bledenharn  and  Mr.  Kenneth  W.  Trahan,  ADP  Center,  are  acknowledged  for 
programming  support. 

Mr.  Donald  R.  Dressier,  was  also  the  OCE  point  of  contact.  The  work  ms 
done  under  the  direction  of  Dr.  N.  Radhakrishnan,  Special  Technical  Assistant, 
ADP  Center.  Mr.  Dressier  and  Dr.  Radhakrishnan  also  contributed  in  the 
definition  of  general  concepts  for  the  development  of  3DSAD.  Mr.  Donald  L. 
Neumann  was  Chief,  ADP  Center. 

Directors  of  WES  during  the  preparation  and  publication  of  this  report 
were  COL  N.  P.  Conover,  CE,  and  COL  T.  C.  Creel,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 
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A  THREE-DIMENSIONAL  STABILITY 
ANALYSIS /DES IGN  PROGRAM  (3DSAD) 

GENERAL  ANALYSIS  MODULE  (CGAM) 

PART  I:  INTRODUCTION 

Overview  of  the  3D  Stability  Program 

1.  The  objective  of  the  Computer-Aided  Structural  Engineering  (CASE) 
Task  Group  on  Three-Dlaenslonal  (3D)  Stability  Analysis  is  to  develop  computer 
programs  that  will  help  design  engineers  perform  stability  computations  for 
general  3D  structures.  To  enable  this,  a  computer  program  called  3DSAD 
(three-dimensional  stability  analysis/design)  is  being  developed  in  a  modular 
fashion.  Initially,  3DSAD  will  have  three  general  modules: 

_a.  General  Geometry  Module.  This  will: 

(1)  Define  geometry  based  on  two-dimensional  (2D)  cross 
sections  extended  in  the  third  dimension,  eight-node 
brick  elements,  or  clusters  of  planar  polygonal 
patches. 

(2)  Perform  centroid,  volume,  and  weight  computations  on 
described  geometry. 

(3)  Employ  interactive  graphics  extensively. 

b.  General  Loads  Module.  This  will  compute  loads  on  general 
3D  structures  based  on  input  of  geometry,  water  levels, 
soil  strata  descriptions,  etc. 

c.  General  Analysis  Module.  This  will  perform  overturning, 
bearing,  and  sliding  computations. 

The  engineer  performing  an  analysis  of  any  3D  structure  will  be  able  to 
interact  directly  with  the  above  modules. 

2.  Besides  the  general  capabilities  that  are  useful  for  any  3D 
structure,  3DSAD  will  also  provide  for  simplified  geometry  and  load  input 
along  with  criteria  check  modules  for  some  specific  structures.  This  latter 
capability  will  permit  interactive  design  of  these  structures.  Examples  of 
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gone  specific  structures  for  which  nodules  will  be  developed  are  dans,  locks 
walls,  powerhouses,  and  pumping  stations. 


3.  A  specific  structure  Input  nodule  requires  less  data  than  that  for  a 
general  structure.  Modules  of  this  type  will  interact  with  the  General 
Geometry  Module  and  the  General  Loads  Module  to  define  the  geometry  and  loads 
internally  in  the  program.  After  analysis,  a  specific  structure  criteria 
check  nodule  will  check  pertinent  values,  change  dimensions  (if  necessary), 
and  cycle  through  the  computation.  A  general  schematic  of  the  3DSAD  program 
is  shown  in  Figure  1. 


Figure  1.  General  schematic  of  3DSAD 


4.  The  3DSAD  program  will  be  developed  in  phases.  During  the  first 
phase,  the  three  general  modules  are  being  developed.  This  approach  will 
enable  the  stability  analysis  of  any  3D  structure  although  the  input  may  be 
more  complicated  than  need  be  for  some  specific  structures.  In  the  subsequent 
phases,  the  special  input  and  criteria  check  nodules  will  be  developed  for 
several  specific  structures. 


The  General  Analysis  Module 


5.  This  report  contains  both  a  user's  guide  and  a  theoretical 
discussion  for  the  General  Analysis  Module.  As  stated  above,  this  nodule 
performs  overturning,  bearing  (base  pressures),  and  sliding  computations  for  a 
structure. 


Input 


6.  Input  to  the  nodule  consists  of: 

a.  The  definition  of  the  area  or  surface  of  investigation 
(generally  the  base  of  the  structure)  by  points, 

line  segments,  and  their  connectivity. 

b.  The  sum  total  of  the  forces  (F^,  Fy,  and  Fz> 

and  moments  (M^,  My,  and  Mz),  excluding  uplift, 

considered  to  compose  the  load  case.  The  forces  and 
moments  are  given  in  consistent  units. 

c.  Uplift  data  consisting  of  uplift  pressure  in  consistent 
units  at  each  point  on  the  base. 

d.  Foundation  data  used  for  computing  the  structure's 
resistance  to  sliding  consisting  of  the  shear  strength 
in  consistent  units  and  the  phi  angle  in  degrees. 


Analysis 


To  do  an  analysis,  the  module: 

a..  Defines  a  plane  of  analysis.  A  modified  base  can  be  used 
which  is  a  planar  approximation  (least  squares  fit)  of 
an  irregular  surface  of  investigation  (original  base) . 

jb.  Using  the  general  flexure  equation,  computes  the  base 
pressures  for  each  point  on  the  base. 

c.  If  any  of  the  pressures  are  negative,  computes  the  new 
effective  base,  modifies  the  uplift  for  category  1  and  2 
load  cases,  and  repeats  the  base  pressure  computation. 
This  process  continues  until  all  the  computed  pressures 
are  positive  (compression)  or  zero. 

d.  Computes  a  shear  friction  factor  of  safety  in  the 
specified  direction. 
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Output 


8.  Output  consists  of: 
a.  A  printout  of: 

(1)  Input  data. 

(2)  Area  properties. 

(3)  Summary  of  forces  and  moments. 

(4)  (X,Y,Z)  coordinates  and  base  pressures 
for  points  defining  the  modified  base. 

(5)  Percent  effective  base  and  shear  friction 
factor  of  safety. 

J>.  Plot  of  outline,  kern,  and  centroid  of  the  base  and 
position  of  the  load  resultant. 


Coordinate  System 


9.  The  coordinate  system  is  shown  in  Figure  2.  Note  that 


Figure  2.  Coordinate  systc 


X  is  to  the  right,  Z  is  up,  and  Y  is  ’’Into”  the  page.  This  is  a  right-handed 
system.  Most  bases  considered  in  the  analysis  will  be  in  the  X-Y  plane 
(horlxontal)  or  nearly  so.  Note  also  that  moments  are  considered 
counterclockwise  positive  (right-hand  screw  rule). 
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PAST  II:  USER'S  GUIDE 


Introduction 

10.  Although  the  General  Analysis  Module  can  be  used  in  conjunction  with 
specific  structures  such  as  dans,  this  input  guide  describes  how  It  can  be 
used  as  a  separate  module.  Data  are  in  the  fora  of  commands  that  the  user 
either  interactively  inputs  or  saves  in  a  data  file  to  be  read  and  executed  by 
the  module  at  one  time.  (See  Appendix  A  for  details  on  creating  and  modifying 
data  files.) 

Running  the  Program 

11.  3DSAD  is  designated  X8100  in  the  Conversationally  Oriented  _Real-Time 
_Program-Generating  JJystem  (CORPS)  library.  Thus,  to  initiate  the  program,  the 
user  types  the  command  on  his  respective  system  that  executes  program  X8100 
from  the  CORPS  library.  The  computer  first  responds  with 


3DSAD  -  CASE 

STRUCTURE  TYPE  OR  GENERAL  MODULE? 


A  ?  response  produces 


A  THREE-DIMENSIONAL  STABILITY  ANALYSIS/DESIGN  PROGRAM 

(3DSAD) 

A  PRODUCT  OF 

Computer  Aided  Structural  Engineering 
(CASE) 


PROGRAM  NO.  713-F3-R0-008 

Computer  Aided  Design  (CAD)  of  STRUCTURAL  STABILITY 

o  o  o  o  o  o 

SR  ?,  HELP,  OR  WHAT  TO  GET  VALID  RESPONSES. 

ER  STOP,  END,  QUIT,  OR  DONE  TO  TERMINATE  PROGRAM. 

ID  SliiuJCTURE  TYPES  ARE:  LOCKS,  DAMS,  POWERHOUSES,  WALLS,  PUMP 
TIONS,  AND  INTAKES. 

ID  GENERAL  MODULE  TYPES  ARE:  GEOMETRY,  LOADS,  AND  ANALYSIS 
UCTURE  TYPE  OR  GENERAL  MODULE  ? 
then  gives 

ANAL 

eneral  Analysis  Module. 

The  next  question  the  program  asks  is 

RESTART  FILE  NAME  OR  CR7 

1"  stands  for  carriage  return.  The  restart  file  saves  all  data  that 
has  thus  far  input  either  from  another  data  file  or  by  typing 
Lvely.  (Note:  utility  commands  such  as  INPUT,  END,  etc.,  are  not 
The  user  giv <s  a  carriage  return  if  he  does  not  want  a  restart  file. 
The  next  question  is 

COMMAND? 

tnda  will  be  discussed  in  detail  in  the  next  section. 


Com  nands 


14.  The  program  uses  commands  (POINTS,  BASE,  CASE,  etc.)  to  build  the 
data.  The  commands  are: 

a.  Data  building: 

POINT  UPLIFT 
CIRCLE  ELLIPSE 
QUADRATIC 
PHI  SHRSTR 
S ANGLE  DENSITY 
BASE  CASE 

b.  Utility: 

INPUT  END 
GO  RETURN 
CLEAR 

This  list  Is  obtained  by  typii«  "?"  at  the  command  level.  Only  the  minimum 
number  of  letters  of  a  command  need  to  be  given.  The  user  can,  however,  type 
the  entire  word  if  he  prefers.  Commands  and  their  accompanying  data  can  be 
put  into  a  data  file  or  typed  interactively  while  the  program  is  running.  All 
lines  in  data  files  created  for  3DSAD  must  be  assigned  line  numbers.  A  file 
can  be  input  at  the  command  level  by  typing  INPUT  and  then  the  file  name. 

This  procedure  is  described  in  detail  in  paragraph  31  below. 

15.  I Each  command  will  now  be  described  in  detail.  In  giving  the  format 
for  the  commands,  actual  letters  to  be  typed  will  be  enclosed  in  quotes  to 
distinguish  them  from  variable  names.  The  quotes  do  not  have  to  be  typed  when 
the  user  issues  the  command.  The  required  letters  are  shown  in  capital 
letters;  the  optional  letters  are  shown  in  lower  letters. 


POINTS 

16.  The  user  first  defines  a  minimum  of  three  points  representing  the 
vertices  of  the  base  using  the  POINTS  command.  Its  format  is 

"POints"  NPT 
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where  MPT  le  Che  number  of  points.  After  this  line,  a  one-  to  four-character 
identification  label  and  the  (X,Y,Z)  coordinates  for  each  point  are  given. 

17.  Interactive  node.  Below  is  an  example  of  the  POINTS  command  done 
interactively  on  a  CDC  system  (all  data  provided  in  the  examples  in  this 
report  have  units  of  feet  for  distance  and  kips  for  force): 


COMMAND  7 
I>POINTS  4 
LABEL,  X,  Y,  Z 
I>1  0  0  268 
t>2  54  0  268 
I>3  54  32  268 
I>4  0  32  268 

18.  Data  file  mode.  The  same  data  are  put  in  a  file  as  shown  below: 


10  POINTS  4 
20  1  0  0  268 
30  2  54  0  268 
40  3  54  32  268 
50  4  0  32  268 

Note  that  these  four  points  could  form  a  horizontal  rectangular  base  at 
elevation  268. 

CIRCLE 

19.  After  the  user  has  defined  some  points  to  work  with,  he  must  then 

define  any  curved  line  segments.  That  is,  line  segments  between  points  are 

assumed  straight  unless  otherwise  specified.  A  circular  arc  may  be  defined  by 

using  the  CIRCLE  command  as  follows: 

"Circle"  N1  N2  R  "Left" 

"Circle"  Nl  N2  R  "Right" 

N1  and  N2  are  two  point  numbers  between  which  a  circular  arc  of  radius  R  is 
being  defined.  "LEFT"  and  "RIGHT"  designate,  respectively,  on  which  side  of 
the  line  segment  the  center  of  the  circle  is  when  progressing  from  Nl  to  N2. 
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If  not  designated,  "LEFT"  is  assumed.  The  following  data  file  results  in  the 
base  shown  in  Figure  3: 


10  POINTS  4 

20  1  0  -5 

30  2  5  0 

40  3  0  5 

50  4  -5  0 

60  CIRC  125 

70  CIRC  2  3  5  LEFT 

80  CIRC  435  RIGHT 

90  CIRC  415 


Figure  3.  CIRCLE  command 


20.  Circles  can  be  defined  only  in  the  X-T  plane. 

ELLIPSE 

21.  The  user  can  define  an  elliptical  line  segment  using  the  ELLIPSE 
command: 

"ELlipse"  N1  N2  A  B  "Left" 

"ELlipse"  N1  N2  A  B  "Right" 

N1  and  N2  are  two  point  labels  between  which  an  elliptical  arc  having 
semimajor  axis  length  A  and  semiminor  axis  length  B  is  drawn.  "LEFT"  and 
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Wf 


"RIGHT"  have  the  ease  meaning  ea  in  the  CIRCLE  coaaand.  The  following  data 


file  reaulta  in  the  base  shown  in  Figure  4: 


10  POINTS  4 
20  1  0  -5 
30  2  10  0 
40  3  0  5 
50  4  '10  0 
60  ELLI  1  2 
70  ELLI  2  3 
80  ELLI  4  3 
90  ELLI  4  1 


10  5 

10  5  LEFT 
10  5  RIGHT 
10  5 


Note  that  the  seaiaajor  and  semioinor  axes  are  always  parallel  to  the 
coordinate  axes. 

22.  As  with  the  CIRCLE  coaaand,  the  ellipse  can  be  defined  only  in  the 
X-Y  plane. 
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QUADRATIC 

23.  Th«  user  uy  need  e  curved  line  segaent  which  Is  not  circular  or 
elliptical.  The  quadratic  line  segaent  Is  provided  for  this  purpose.  The 
caaaand  forest  Is 

"Quadratic"  Nl,  N2,  XQQ,  YQQ,  ZQQ 

N1  end  N2  ere  the  point  labels  of  the  beginning  and  end  of  the  quadratic  line 
segaent  and  (XQQ,  YQQ,  ZQQ)  Is  an  Interpolation  point  (Figure  5)  that  the 


(XQQ,  YQQ,  ZQQ) 


Figure  5.  Quadratic  plotting 


curve  aust  pass  through.  The  following  data  file  results  In  the  base  shown  In 


Figure  6: 

10  POINTS  4 
20  1  0  -5 
30  2  5  0 
40  3  0  5 
50  4  -5  0 
60  QUAD  124-3 
70  QUAD  2334 
80  QUAD  34-43 
90  QUAD  41-3-4 
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24.  The  user  next  describes  the  order  in  which  the  points  and  curved 
line  segments  are  connected  by  giving  the  BASE  cooaand.  The  format  for  this 
command  is 

"Base"  NHOLE 

where  NHOLE  is  the  number  of  holes  in  the  base.  After  the  BASE  command  is 
given  in  the  interactive  mode,  the  connectivity  data  are  requested  for  the 
outer  boundary  and  each  hole: 

CONNECTIVITY  DATA? 

The  connectivity  data  (see  Pigure  7)  are  described  by  first  giving  the  number 


Figure  7.  Sample  base 


of  points  forming  the  outer  boundary  followed  by  the  point  labels  in 
counterclockwise  order  for  the  outer  boundary  and  then  the  number  of  points 
forming  the  hole  followed  by  the  point  labels  in  clockwise  order.  The 
interactive  question  and  answer  sequence  used  to  generate  the  base  in  Figure  7 


COMMAND  ? 

DPOINTS  8 
LABEL  X,  Y,  Z 
I>1  a  e  268 
I>2  54  a  268 
I>3  54  32  268 
I>4  d  32  268 
I?5  2d  Id  268 
I>6  3d  Id  268 
I>?  3d  2d  268 
I>8  2d  2d  268 
COMMAND  ? 

DBASE  1 

OUTER  BOUNDARY  DATA 
CONNECTIUITY  DATA  ? 

1)4  1  2  3  4 

DATA  FOR  HOLE  1 

CONNECTIUITY  DATA  ? 

I>4  8  7  6  5 


If  these  data  were  placed  into  a  data  file,  the  file  would  appear  as: 


Id  POINTS  8 
30  1  d  d  268 
3d  2  54  d  268 
4d  3  54  32  268 
5d  4  d  32  268 
60  5  2d  id  268 
70  6  3d  id  268 
80  7  3d  20  268 
90  8  2d  2d  268 
100  BASE  1 
lie  4  1.2  3  4 
12d  4  8  7  6  5 

UPLIFT 


25.  When  one  or  more  of  the  base  pressures  are  negative  on  the  first 


Iteration,  a  new  effective  base  and,  possibly,  a  new  value  for  uplift  (see 


Part  III)  oust  be  coaputed.  Thus,  the  uplift  data  oust  be  specified 


separately  froa  the  other  forces  and  aoaents  of  the  load  case.  These  data  are 


provided  by  giving  the  uplift  pressure  at  each  of  the  defined  points  on  the 


base.  The  format  for  the  UPLIFT  command  is 


"Uplift"  NPT 

where  NPT  is  the  number  of  points  for  which  uplift  is  now  to  be  specified. 
After  this  line,  the  identification  label  and  corresponding  value  of  uplift 
pressure  are  given  for  each  point.  Below  Is  an  exaaple  of  the  UPLIFT  command: 

10  POINTS  4 
20  1  0  0  268 
30  2  54  0  268 
40  3  54  32  268 
50  4  0  32  268 
60  UPLIFT  4 
70  1  2.000 
80  2  0.125 
90  3  0.125 
100  4  2.000 

Zero  is  assumed  for  uplift  for  all  points  not  specified. 

PHI 

26.  The  resistance  to  sliding  is  computed  from  the  formula 

R  -  N  tan*  +  sA  (1) 

where  R  is  the  resistance,  N  is  the  normal  component  of  the  force  applied  to 
the  base,  4  (phi)  is  the  angle  of  Internal  friction  of  the  foundation  material 
in  degrees,  s  is  the  shear  strength  or  cohesion  of  the  foundation  material, 
and  A  is  the  effective  area.  The  PHI  command  is  used  to  specify  the  4 
angle.  Its  format  is 

"PHI"  ANGLE 

where  ANGLE  is  an  angle  specified  in  degrees. 

SHRSTR 

27.  THE  SHRSTR  command  is  used  to  specify  the  shear  strength  in 
consistent  units.  Its  format  is 
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"SHrstr"  S 

where  S  Is  the  shear  strength. 

SANGLE 

28.  The  S ANGLE  command  allows  the  user  to  specify  in  what  direction  the 
shear  friction  factor  of  safety,  as  given  in  Engineer  Technical  Letter  (ETL) 
1110-2-184,*  is  to  be  computed.  Its  format  is: 

"SAngle"  "Default" 

"SAngle"  SANG 

SANG  is  the  angle  in  the  horizontal  plane  measured  from  the  global  X  axis  for 
which  sliding  stability  is  to  be  evaluated.  The  default  angle  is  the  angle 
formed  by  the  horizontal  component  of  the  resultant  load  and  the  global  X 
axis.  Note,  however,  that  sliding  is  always  computed  tangent  to  the  plane  of 
analysis. 

DENSITY 

29.  The  default  value  for  the  density  of  water  is  assumed  to  be  0.0625 
kipa/ft^  (0.001  kg/m^).  The  user  can  change  this  value  to  his  set  of 
consistent  units  by  using  the  DENSITY  command.  Its  format  is 

"Density"  DENS 

where  DENS  is  the  specified  density  of  water.  Note  that  this  command  would  be 
used  if,  for  instance,  a  metric  system  of  units  were  being  used. 

CASE 

30.  After  the  user  has  specified  points  and  any  curves,  defined  the 
connectivity  for  the  base,  specified  any  uplift  values,  and  defined  PHI, 
SHRSTR,  and  DENSITY  when  necessary,  he  must  give  the  CASE  command  to  define 


*  This  and  other  engineer  standards  used  in  development  of  the  General 
Analysis  Module  are  listed  in  the  Bibliography. 
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Che  other  applied  forces  and  moments.  ics  format  is 


‘‘Case”  CNAME ,  LDTYPE 

where  CNAME  Is  a  one-  Co  four-character  identifier  and  LDTYPE  is  Che  load 
cype.  LDTYPE  may  have  Chree  values: 

a.  1  -  long  cera. 

b.  2  -  shore  cerm. 

3  -  instantaneous. 

Specifying  an  insCantaneous  load  Cype  (LDTYPE-3)  prevents  Che  uplift  from 
being  modified  when  Che  effective  base  is  modified.  After  the  CASE  command 
is  given,  Che  next  data  required  (either  interactively  or  from  a  data  file) 
are  the  three  force  components  (Fx,  Fy,  Fz)  and  moments  (Mx,  My,  M^)  about  the 
global  coordinate  system  excluding  the  contribution  from  uplift.  An  example 
data  file  is 


ie  POINTS  4 
20  1  0  0  268 
30  2  54  0  268 
40  3  54  32  268 
50  4  0  32  268 
60  UPLIFT  4 
70  1  2.000 
80  2  0.125 
90  3  0.125 
100  4  2.000 
110  PHI  19.6 
120  SHRSTR  .22 
130  BASE 
140  4  1  2  3  4 
150  CASE  NORN  1 

160  1100  0  -5800  -93000  444000  0 


After  the  date  for  the  CASE  command  are  given,  the  program  performs  the 
computations  and  begins  output.  A  detailed  discussion  of  the  output  begins  in 
paragraph  36. 

INPUT 

31.  This  command  allows  the  user  to  input  or  read  into  memory  a 
permanent  data  file  saved  on  disc.  Its  format  is 


"INPut"  FLNM1 
"INPut"  FLNM1  "P" 

where  FLIM1  is  a  1-  to  20-character  file  description.  If  the  P  is  also  typed, 
a  detailed  printout  of  the  input  file  is  printed  on  the  terminal  as  if  it  had 
been  done  interactively. 

END 

32.  This  command  is  given  to  terminate  running  of  the  program.  Its 
format  is 


"ENd" 


GO 

33.  This  command  is  used  when  the  program  is  being  used  for  a  specific 
structure  such  as  a  dam  or  lock.  Giving  the  GO  command  automatically  causes 
the  program  to  return  to  do  the  next  load  case.  The  format  for  this  command 
is 

"GO" 

RETURN 

34.  This  command  is  used  to  return  to  the  question 


STRUCTURE  TYPE  OR  GENERAL  MODULE? 


so  that  the  user  can  select  another  nodule.  The  format  for  this  command  Is 


"REturn" 

CLEAR 

35.  This  command  Is  used  to  clear  all  defialcioa  of  geometry  from  memory 
so  as  to  begin  a  new  problem.  Its  format  is 

"CLear" 


Output 


36.  Each  time  a  CASE  command  is  given,  an  analysis  is  performed  and  the 
output  immediately  follows. 

Printout  of  input 

37.  The  user  is  first  asked: 

DO  YOU  WISH  TO  SEE  A  SUMMARY  OP  THE  ANALYSIS  INPUT  DATA? 

If  YES  Is  typed,  the  following  table  is  printed: 

LOAD  CASE  NORM  CATEGORY  *  1 


INPUT  FORCES  AND  MOMENTS  EXCLUDING  UPLIFT 


FX  * 

1100.000 

FY  *• 

0. 

FZ  * 

-5800.000 

MX  ■ 

-93000.0 

MY  « 

444000.0 

MZ  * 

0. 

INPUT 

COORDINATES  AND  UPLIFT 

VALUES 

NAME 

X 

Y 

Z 

UPLIFT 

1 

0. 

0. 

268.000 

2.000 

2 

54.000 

0. 

268.000 

0.125 

3 

54.000 

32.000 

268.000 

0.125 

4 

0. 

32.000 

268.000 

2.000 

PHI  - 

19.60  8HR8TR  > 

0.22 

SANGLE  * 

DEFAULT 

j 

I 


DENSITY  * 


0.0625 


Final  kern  plot 

38.  The  user  in  nose  cases  is  next  asked  the  question: 

DO  YOU  WISH  TO  SEE  THE  FINAL  KEEN  PLOT? 

IF  YES  is  given,  the  following  are  plotted: 

a.  Outline  of  the  modified  base  (original  base  projected 
onto  the  plane  of  investigation)  with  solid  lines. 

t>.  Kern  with  dashed  lines. 

c^.  Circled  +  marking  the  centroid  of  the  base, 
d^.  +  marking  the  position  of  the  load  resultant. 

_e.  Values  of  eccentricity. 

f_.  Distance  of  load  resultant  to  kern. 

Figure  8  shows  the  final  kern  plot  for  the  data  file  given  in  paragraph  30 
When  this  plot  is  completed,  the  bell  will  ring  and  the  program  will  pause 


Figure  8.  Final  kern  plot 


:s 


>  continue,  the  user  hits  the  carriage  ret^*n.  The  final  results 
>llows : 


LOAD  CASE  NORM  CATEGORY  =  1 


FINAL  BASE  AREA  PROPERTIES 


1728.000 

IXP  =  147456.0 

IYP  « 

419904.0 

27.000 

YBAR  =  16.000 

ZBAR  = 

268.000 

0. 

ZNANG  =  0. 

PXANG  = 

0. 

SUMMARY 

OF  FINAL  FORCES  AND 

MOMENTS 

—  INPUT— 

1100.000 

FY  =  0. 

FZ  = 

-5800.000 

-93000.0 

MY  =  444000.0 

MZ  =• 

0. 

—COMPUTED  UPLIFT  — 

0. 

FY  =  0. 

FZ  = 

1836.000 

29376.0 

MY  =  -34992.0 

MZ  = 

0. 

— TOTAL— 

1100.000 

FY  =  0. 

FZ  = 

-3964.000 

-63424 . 0 

MY  =  409008.0 

MZ  = 

0. 

FINAL  IN  PLANE  COORDINATES  AND  BASE  PRESSURES 

X 

Y 

Z 

PRESSURE 

0. 

0. 

268.000 

1.811 

54.000 

0. 

268.000 

2.734 

54.000 

32.000 

268.000 

2.777 

0. 

32.000 

268.000 

1.854 

FINAL  SHEAR  FRICTION  FACTOR  OF 

SAFETY 

19.60  SHRSTR  =  0.22 

SANGLE  * 

0. 

FACTOR  OF  SAFETY  =  1.63 


EFFECTIVE  BASE  *  100. OX 


le  output  data  are  described  In  the  following  tabulation: 
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Variable 


CATEGORY 

AREA 

IXP 

IYP 

XBAR 

YBAR 

ZBAR 

XYANG 

ZNANG 

PXANG 

FX 

FY 

FZ 

MX 

MY 

MZ 

NAME 

(X.Y.Z) 

PRESSURE 

FACTOR  OF 
SAFETY 


_ Description _ 

1  ~  long  term,  2  -  short  term,  3  -  instantaneous 

Area  of  the  base  in  the  plane  of  analysis 

X  principal  moment  of  area  in  the  local 
coordinate  system 

Y  principal  moment  of  area  in  the  local 
coordinate  system 

X  centroid  in  global  coordinates 

Y  centroid  in  global  coordinates 

Z  centroid  in  global  coordinates 

Angle  in  the  X-Y  plane  that  the 
direction  of  the  plane  of  analysis  makes 
with  the  global  X-axis  (see  Figure  9).  This 
corresponds  to  a  rotation  about  the  global  Z-axis 
to  form  the  prime  set  of  axes 

Angle  between  the  global  Z  axis  and 
the  upward  normal  to  the  plane  of 
analysis  (see  Figure  9). 

This  corresponds  to  a  rotation  about  the  Y'-axis 
to  form  the  double  prime  coordinate  system 

Angle  between  local  X  axis  and  the 
principal  X  axis  in  the  local  coordinate 
system  (see  the  angle  THETA  in  Figure  15) 

Force  in  the  global  X  direction 

Force  in  the  global  Y  direction 

Force  in  the  global  Z  direction 

Moment  about  the  global  X  axis 

Moment  about  the  global  Y  axis 

Moment  about  the  global  Z  axis 

Point  label 

In  plane  coordinates  of  point 
Computed  bearing  or  base  pressure 
Shear  friction  factor  of  safety 
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Pp'-allel  to  downwardmost 
direction 


Figure  9.  Angles  of  rotation 


Intermediate  kern  plot 

41.  If  any  of  the  base  pressures  are  negative  and  the  program  has  to 
iterate  to  find  a  new  effective  base,  the  user  first  has  the  option  of  viewing 
a  kern  plot  before  the  iteration  begins.  The  question  asked  is 
DO  YOU  WANT  TO  SEE  THE  INTERMEDIATE  KERN  PLOT? 

An  example  of  the  intermediate  kern  plot  is  given  in  paragraph  44. 
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Intermediate  results 


42.  When  the  intermediate  kern  plot  is  complete,  the  bell  will  ring  and 
the  program  will  pause.  A  carriage  return  will  cause  the  program  to  continue 
and  print  the  Intermediate  results  (see  the  example  in  paragraph  44).  If  all 
the  intermediate  base  pressures  are  zero  or  positive,  the  final  results  are 
the  same  as  the  intermediate  ones  and  thus  are  the  only  ones  printed. 

Note  also  that  the  load  resultant  is  inside  the  kern  when  the  base  pressures 
are  all  positive. 


Additional  Sample  Problems 


43.  To  help  illustrate  the  features  of  this  module,  some  additional 
sample  problems  will  now  be  presented. 

Negative  pressures 

44.  The  following  data  file  is  a  simple  example  of  how  the  iterative 

process  works  when  one  or  more  of  the  initial  base  pressures  are  negative: 

10  POINTS  4 
20  1  0  0 
30  2  10  0 
40  3  10  10 
50  4  0  10 
60  BASE 
70  4  1  2  3  4 
80  CASE  TEST  1 
90  0  0  -1  -7  7  0 


The  intermediate  kern  plot  is  shown  in  Figure  10. 
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Figure  10.  Intermediate  kern  plot, 
negative  pressures  example 


The  intermediate  results  are: 

LOAD  CASE  TEST  CATEGORY  -  1 


COMPUTED  BASE  AREA  PROPERTIES 


AREA  =  100.000  IXP 

XBAR  =  5.000  YBAR 

XYAMG  *  0.  ZNANG 


833.3  IYP  * 
5.000  ZBAR  - 
0.  PXANG  = 


FX  = 
MX  = 


FX  = 
MX  * 


FX  * 
MX  * 


SUMMARY  OF  FORCES  AND  MOMENTS 
—  INPUT— 

0.  FY  =  0.  FZ  = 

-7.0  MY  =  7.0  HZ  * 

—COMPUTED  UPLIFT— 

0.  FY  *  0.  FZ  = 

0.  MY  =  0.  MZ  - 

—TOTAL— 

0.  FY  =  0.  FZ  * 

-7.0  MY  =»  7.0  HZ  » 

(Continued) 


833.3 

0. 

0. 


-1.000 

0. 


0. 

0. 


-1.000 

0. 
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COMPUTED  IN  PLANE  COORDINATES  AND  BASE  PRESSURES 


NAME 

X 

Y 

Z 

PRESSURE 

1 

0. 

0. 

0. 

-0.014 

•» 

a- 

10.000 

0. 

0. 

0.010 

3 

10.000 

10.000 

0. 

0.034 

4 

0. 

10.000 

0. 

0.010 

COMPUTED  SHEAR 

FRICTION  FACTOR 

OF  SAFETY 

PHI  * 

0.  SHRSTR 

=•  0. 

SANGLE  * 

0. 

FACTOR  OF  SAFETY 

*  9999.99 

The  final  kern  plot  Is  given  In  Figure  11. 


Figure  11.  Final  kern  plot 
negative  pressures  example 


The  results  after  Iteration  are: 


LOAD  CASE  TEST  CATEGORY  -  1 


FINAL  BASE  AREA  PROPERTIES 


AREA  = 

68.846 

IXP  - 

271.6 

IYP  - 

XBAR  * 

6.072 

YBAR  * 

6.072 

ZBAR  - 

XYANG  * 

0. 

ZNANG  - 

0. 

PXANG  - 

SUMMARY 

OF  FINAL  FORCES  AND 

MOMENTS 

-INPUT-- 

FX  * 

0. 

FY  = 

0. 

FZ  - 

MX  = 

-7.0 

MY  =• 

7.0 

HZ  * 

—  COMPUTED  UPLIFT— 

FX  = 

0. 

FY  -- 

0. 

FZ  * 

MX  =• 

0. 

MY  =• 

0. 

MZ  * 

-TOTAL-- 

FX  = 

0. 

FY  * 

0. 

FZ  = 

MX  * 

-7.0 

MY  " 

7.0 

MZ  => 

FINAL  IN  PLANE 

COORDINATES  AND 

BASE  PRESSURES 

NAME 

X 

Y 

Z 

2-1 

7 .894 

0. 

0. 

2 

10.000 

0. 

0. 

3 

10.000 

10.000 

0. 

4 

0. 

10.000 

0. 

4-1 

0. 

7.894 

0. 

FINAL  SHEAR  FRICTION  FACTOR  OF  SAFETY 
PHI  *  0.  SHR8TR  -  0.  SANGLE  - 

FACTOR  OF  SAFETY  ■  9999.99 


671.6 

0. 

-45.000 


-1.000 

0. 


0. 

0. 


1.000 

0. 


PRESSURE 

-0.000 

0.007 

0.041 

0.007 

0.000 


0. 


EFFECTIVE  BASE 


68.8% 


Honhorlzontal  nonplanar  bate 

45.  The  Input  date  are  aa  follows: 

10  P01N  4 
20  1  0  0  -1 
30  2  10  0  0 
40  3  10  10  2 
50  4  0  10  0 
60  BASE 
70  4  1  2  3  4 
80  CASE  Cl  1 
90  0  0  -1  -5.5  5.5  0 

Figure  12  shows  the  final  kern  plot . 


Figure  12.  Final  kern  plot,  nonhorizontal 
nonplanar  base  example 


Final  results  are 


LOAD  CASE  Cl  CATEGORY  -  1 


FINAL  BASE  AREA  PROPERTIES 


AREA  - 

102.225 

IXP  *  890.3 

IYP  * 

851.9 

XBAR  ■ 

5.012 

YBAR  *  5.012 

ZBAR  = 

0.25* 

XYANG  » 

-135.000 

ZNANG  =  11.977 

PXANG  = 

-45.000 

SUNNARY 

OF  FINAL  FORCES  AND 

N0NENT8 

— INPUT — 

FX  » 

0. 

FY  =  0. 

FZ  * 

-1.000 

NX  « 

-5.5 

NY  *  5.5 

NZ  * 

0. 

—  CONPUTED  UPLIFT— 

FX  * 

0. 

FY  =  0. 

FZ  * 

0. 

NX  * 

0. 

NY  *  0. 

NZ  - 

0. 

—TOTAL— 

FX  * 

0. 

FY  *  0. 

FZ  = 

-1.000 

NX  - 

-5.5 

NY  =■  5.5 

NZ  * 

0. 

FINAL  IN  PLANE  COORDINATES  AND  BASE  PRESSURES 

NANE 

X 

Y 

Z 

PRESSURE 

1 

0.036 

0.036 

-1.239 

0.004 

2 

9.964 

-0.036 

0.239 

0.010 

3 

10.036 

10.036 

1.761 

0.015 

4 

-0.036 

9.964 

0,239 

0.010 

FINAL  SHEAR  FRICTION  FACTOR  OF 

SAFETY 

PHI  = 

0.  SHR8TR  =  0. 

SANGLE  = 

-135.00 

FACTOR  OF  SAFETY  ■  0. 

EFFECTIVE  BASE  -  100.01 


Powerhouse  (uulyiii 

46.  The  following  example  is  an  analysis  of  a  powerhouse  provided  by  Mr 
Charles  W.  Kling,  Mobile  District.  The  data  file  for  the  analysis  appears  as 
follows : 


10  P0  22 
20  D  80  113 
30  C  0  113 
40  B  0  0 
SO  A  80  0 
60  1  72.5  113 
70  2  72.5  83.4 
80  3  71  73.5 
90  4  65.7  63 
100  5  59.2  53 
110  6  52.6  39.9 
120  7  34.2  38 
130  8  21.9  S3 
140  9  13.2  63 
150  10  9  73.5 
160  11  7.5  83.4 
170  12  7.5  113 
180  13  36.5  113 
190  14  36.5  73.5 
200  15  39.1  63 
210  16  42.6  63 
220  17  43.5  73.5 
230  18  43.5  113 
240  BASE  1 
250  4  A  B  C  D 

260  18  1  2  3  4  5  6  7  8  9  10  11  12 
270  13  14  15  16  17  18 

**80  CASE  1ST  2 

290  0  0  -35562  -2126942  1473689  0 


The  final  kern  plot  is  shown  in  Figure  13.  The  output  appears  as  follows: 

LOAD  CASE  1ST  CATEGORY  =  2 


FINAL  BASE  AREA  PROPERTIES 


AREA  = 
XBAR  - 
XYANG  * 


5374.850 

39.934 

0. 


IXP  = 
YBAR  - 
ZNAN6  * 


4684220.6 

40.255 

0. 


IYP  * 
ZBAR  - 
PXAN6  » 


3661498. 

0. 

-0.13 


(Continued) 
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i 


Figure  13.  Final  kern  plot, 
powerhouse  analysis  example 
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SUMMARY  OF  FINAL  FORCES  AND  MOMENTS 


-INPUT- 


FX  - 
MX  - 

0. 

-2124942.0 

FY  *  0. 

MY  =•  1473489.0 

F2  = 
MZ  - 

-35542.000 

0. 

—  COMPUTED  UPLIFT— 

FX  ~ 
MX  * 

0. 

0. 

FY  =  0. 

MY  =■  0. 

FZ  -- 
MZ  = 

0. 

0. 

— TOTAL— 

FX  * 
MX  - 

0. 

-2124942.0 

FY  =  0. 

MY  -  1473489.0 

FZ 

MZ  - 

-35542.000 

0. 

FINAL  IN  PLANE 

COORDINATES  AND 

BASE  PRESSURES 

NAME 

X 

Y 

Z 

PRESSURE 

A 

80.000 

0. 

0. 

1.223 

B 

0. 

0. 

0. 

0.040 

C 

0. 

113.000 

0. 

14.835 

D 

80.000 

113.000 

0. 

17.997 

1 

72.500 

113.000 

0. 

17.888 

2 

72.500 

83.400 

0. 

13.494 

3 

71.000 

73.500 

0. 

12.003 

4 

45.700 

43.000 

0. 

10.347 

5 

59,200 

53.000 

0. 

S.7S8 

4 

52.400 

39.900 

0. 

4.748 

7 

34.200 

38.000 

0. 

4.198 

3 

21.900 

53.000 

0. 

8.244 

9 

13.200 

43.000 

0. 

9.404 

10 

9.000 

73.500 

0. 

11.102 

11 

7.500 

83.400 

0. 

12.550 

12 

7.500 

113.000 

0, 

14.944 

13 

34.500 

113.000 

0. 

17.345 

14 

34.500 

73.500 

0. 

11.502 

15 

39.100 

43.000 

0. 

9.981 

14 

42.400 

43.000 

0. 

10.032 

17 

43.500 

73.500 

0. 

11.403 

18 

43.500 

113.000 

0. 

17.447 

FINAL  SHEAR  FRICTION  FACTOR  OF  SAFETY 
PHI  =  0.  SHRSTR  *  0.  SANOLE  *  0. 

FACTOR  OF  SAFETY  *  9999.99 


EFFECTIVE  BASE  *  100. OX 


PART  III:  THEORY  AND  PROCEDURE 


Introduction 


47.  This  part  describes  how  some  of  the  more  significant  output  data  are 
computed.  The  description  will  involve  a  theoretical  discussion  in  some 
cases,  such  as  that  of  base  pressure  computation.  Other  times,  such  as  for 
the  description  on  how  uplift  is  computed,  the  discussion  will  take  the  form 
of  a  description  of  procedure. 

Assumptions 


48.  The  following  assumptions  are  made  in  thl»  analysis: 

j*.  The  effects  due  to  the  elastic  properties  of  both 
the  structure  and  foundation  have  been  ignored; 
rigid  body  motion  assumed;  and  limit  equilibrium 
methods  utilized. 

_b.  The  soil  and  water  loads  are  computed  separately 
from  that  of  the  structure.  Thus,  no  soil- 
structure  interaction  is  considered. 

c_.  The  analysis  is  performed  on  a  planar  approximation 
of  the  base. 

d.  The  analysis  is  thus  performed  using  only  the 
force  components  (F^,  F  ,  ?z)  and  moment 
components  (Mx,  M^,  Mz)  and  the  planar 
approximation  of  the  base. 

je.  The  bearing  or  base  pressures  thus  vary  linearly  over 
the  base  (general  flexure  formula  for  biaxial  bending). 

_f.  The  shear  friction  formula  is  used  for  computing 
the  factor  of  safety  against  sliding. 
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Sign  Convention 

49.  A  strictly  right-handed  coordinate  system  for  coordinates,  forces, 
and  moments  is  used.  The  weight  of  the  structure,  therefore,  is  considered  a 
negative  number.  Compressive  stresses  or  pressures  are  considered  positive, 
while  tensile  stresses  are  negative. 

Base  Pressures 


Formula 

50.  Consider  the  planar  horizontal  base  with  the  origin  at  the  centroid 
as  shown  in  Figure  14.  The  base  pressures  are  computed  from  the  general 

y 


.X 


Figure  14.  Horizontal  base 


flexure  equation 


if  the  principal  axes  of  the  base  are  also  the  X  and  Y  axes.  Here,  A  is  the 
area  of  the  base,  I  is  the  X  moment  of  area,  and  I  is  the  Y  moment  of  area. 

*  y 
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,et  us  start  with  the  assumption  that  pressure  varies  linearly  over 
Thus, 

p  -  a  +  bX  +  cY  (3) 

,  and  c  are  constants  to  be  evaluated. 

First,  integrate  both  sides  of  Equation  3  over  the  area  of  the  base. 
/Apda  »  a  da  +  b  J^Xda  +  c  J^Yda 


ling  Z  forces, 


/Apda  +  Fz  -  0 

/^Xda  *  0 
/A Yda  -  0 

X  and  Y  centroid  are  both  zero.  Therefore, 
-Fz  -  aA  +  b(0)  +  c(0) 


Second,  multiply  both  sides  of  Equation  3  by  X  and  integrate  over 


/AXpda  -  a/AXda  +  b/AX2da  +  c^XYda 


using  Y  momenta, 


-  /AXpda  +  My  -  0 
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By  definition 


*» * ■ 


» 


S '  V2** 


and  since  the  principal  axes  coincide  with  the  coordinate  axes. 


/AXYda  -  0 


Thus, 


My  -  a(0)  +  bly  +  c(0) 


54.  Last,  multiply  both  sides  of  Equation  3  by  Y  and  Integrate  over  the 


area: 


/AYpda  -  a  J^Yda  +  b  ^XYda  +  c  ■rAY2da 

Summing  X  moments. 


By  definition, 


/AYpda  +  Mx  -  0 


Ix  "  'A*2da 

Collecting  terms, 

-Mx  -  a(0)  +  b(0)  +  c  Ix 
M 

x 
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Combining  Equations  3,  4,  5,  and  6  gives  the  desired  result: 

F 

z 

p  - 

General  formulation 

55.  In  general,  however,  the  coordinate  system  will  neither  have  Its 
origin  at  the  centroid  of  the  base  nor  align  with  its  principal  axes  (see 
Figure  15).  The  X'-Y'  coordinate  system  has  its  origin  at  the  centroid  of  the 
base  and  aligns  with  the  X-Y  coordinate  system.  The  X"-Y"  coordinate  system 


'Figure  15.  General  base 


is  aligned  with  the  principal  axes  of  the  plane  and  is  achieved  by  rotating 
the  X'-Y’  axes  an  angle  ©.  Thus,  from  Equation  2, 


F  /M"\  /M"\ 

,m'*wW 

To  put  this  equation  in  terns  of  the  X-Y  system,  both  a  translation  and  a 
rotation  are  required.  First, 


X'  -  X  -  xc 
Y'  -  Y  -  Y. 


Next, 


[X“  Y"]  -  [X'  Y'l 


080  -sinl 
in©  cosi 


Applying  this  translation  and  rotation  to  the  moaent  gives 

<  -  Mx  -  F*Yc 


m;  -  My  +  FbXc 


t"  mr  -  Fm’  r« 

'  ]  ML* 


cos©  -sin© 
sin ©  cos© 


Applying  this  translation  and  rotation  to  the  moment  of  areas  gives 

Xx"  ;AY’2da 

-  ;A  (Y  -  Yc)2da 

-  Ix  -  2YC(YCA)  +  Y£2(A) 

"  lx-YcA 
Xxy  *  -^X'Y'da 

“  /A  (*  -  *C)(Y  -  Yc)da 

"  xXy  "  Xc<YcA>  ‘  Yc<XcA>  +  XcYcA 
“  Xxy  ”  XcYcA 
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x; 


Equation  7  can  now  be  written 


/Ax'2da 


-  /X(X  -  Xc)2da 

-  ly  -  2Xc(XcA)  4  X*A 


Iy  -  X‘A 


f  a-  [■ 


COSO  Bind 
•alnd  coed 


x  ”Ixy  co#®  -•ln* 

f  Bind  coed 


LXxy 


:] 


p  ■  - 


l£  +  [X“XC  Y-YJ 


coed  -Bind 
sind  cosd 


To  get  the  font 


p  ■  a 


+  bX  +  cY 
F 


a  -  -  x"  ♦  [Xc  Yc]  |  coed  -Bind 


coed  -elnd 
8lnd  coed 


M" 

-JL 

I" 

y 

M" 

x 

I" 

X 


(2ble) 


M” 

JL 

I" 

y 


M" 
x 

I" 

X  -I 


M"  M" 

y  x 

b  »  jn  coed  +  Y"  Bio® 

y  x 


(8) 


M" 


M” 


c  **  8lnd  -  -JT7 


coed 


Kern  Computation 


Definition 


56.  Given  the  base  in  Figure  16,  let  P  be  the  point  at  which  the  normal 
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Figure  16.  Sample  base 


force  F  is  applied  to  the  base  to  produce  the  moments  M  and  M •  Thus, 
z  *  y 

Mx  -  Fz  (Yc  +  L  sin  o)  (9) 

My  -  -Fz  <Xc  +  L  cos  a) 

The  kern  is  defined  as  that  area  where,  if  the  point  P  falls  inside  it,  all  of 
the  computed  base  pressures  are  positive.  This  means  that  if  P  falls  outside 
the  kern,  some  points  on  the  base  will  have  a  negative  value. 


Derivation 


57.  Rather  than  start  with  Equation  8,  let  ua  derive  a  new  set  of  a,  b, 
and  c  for 

p  -  a  +  b(X  -  Xc>  +  c(Y  -  Yc)  (10) 

First,  define  the  constants 

\  •  'am* 

\  -  /AU. 

<U) 

V  " 

V  V2d* 


Integrating  both  sides  of  Equation  10  as  before  gives 

/A pda  -  a  /Ada  +  b  /A(X  -  Xc)da  +  c  ^(Y  -  Yc)da 
-Fz  -  aA  +  b(Ax  -  XCA)  +  c(Ay  -  Y£A) 

■  aA 


F 

/X pda  -  ~  a“  j[Xda  +  b^A(X  -  Xc)Xda  -  Y£)Xda 


My  "  'Fr  T  +  b(Axx  “  W  +  C<V  -  V*> 

A  A^ 

H,  -  -F,  J  +  KA„  -  +  c(Ajcy  -  Vjr>  (12) 
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/Ap¥da  •  —  — —  ^Yda  +  b  -^(X  -  Xc)Yda  c  ^(Y  -  Yc)Yda 


A 


-Mx  -  A  +  ***y  -  XcAy>  +  c(Ayy  '  Vy> 


A  A 


-M  -  -F  _Z  +  b(A  -  )  +  c(A  -  ) 

*  z  a  xy  A  w  A 


yy  A 


Combining  Equations  9-13  gives 


-F_ 


- 

p-  — 

r  1 

Xc 

-F  L 
z 

sin  a 

F 

z 

Yc 

cob  a 

A 

i _ 

A  A 


^xx  ~  a  A*y  “ 


A  A  A2 

A  _  A  _  JL 

xy  A  \y 


Collecting  terms. 


A2  A  A  ' 

r  ■ 

*  -  x  *  -  x  y 

*xx  a  ^y  A 

:  b 

i 

m  -f  L 

Z 

sin  a 

? 

cos  a 

A  A  A2 

i 

““  * 

a  _  x  y  a _ * 

L^y  A  yy  a 

I  C 

or 


--FZL 


[sin  a  | 
cos  aj 


-  T  FZL 


Putting  these  results  into  Equation  10  gives 


P  -  -  ^  f1  +  CX-Xc  Y-Yc]  l  1 


\ 

sin  a 

cos  a 

) 
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(13) 

i  r  i 


wmmm 


m 


►  * 


58.  The  kern  is  computed  using  the  following  two  assumptions: 

ju  The  force  Fz  applied  at  any  point  on  the  kern  will 

cause  the  pressures  to  be  zero  at  some  point  on  the  base. 

b.  Assuming  the  base  is  approximated  by  a  series  of  straight 

lines,  the  point  where  pressure  first  becomes  zero  will  always 
be  a  vertex  of  the  polygon  forming  the  base. 

Thus,  along  any  ray  forming  an  angle  a  with  the  horizontal,  find  the  length  L 

such  that  all  pressures  at  the  vertices  of  the  base  are  zero  or  positive.  The 

value  of  L  making  the  pressure  zero  at  the  “ith"  vertex  from  Equation  14  is 


(15) 


AiVxc  W  =  1  f8inot] 

[_cos  aj 

Check  all  vertices  and  choose  the  smallest  positive  L^. 


Lkern  *  Z*  111111  (Positive  L±> 


This  process  is  repeated  until  the  kern  is  clearly  defined. 

Base  of  Analysis 


59.  All  of  the  above  derivations  start  with  the  assumption  of  a 
horizontal  base.  Actually,  the  user  specifies  an  area  or  surface  of 
investigation  which  may  be  neither  horizontal  nor  planar.  What  the  program 
does,  therefore,  is  to  compute  a  planar  base  of  analysis  from  the  base  input 
by  the  user. 

Plane  of  analysis 

60.  The  plane  of  analysis  is  computed  by  doing  a  least  squares  fit  on 
the  N  data  points  describing  the  base.  The  criterion  used  is  to  minimize  the 
sum  of  the  vertical  distance  squared  from  the  N  points  to  the  plane  of 
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analysis  to  be  determined.  The  equation  of  the  plane  is 


Z  -  a  +  bX  +  cY  (16) 

where  a,  b,  and  c  are  constants  to  be  computed*  The  least  squares  expression 
is 


i?1  (a  +  bX±  +  cY1  -  Zt)2  -  nin 

This  gives  the  following  systea  of  simultaneous  linear  equations  to  solve  for 
a,  b,  and  c: 


N  N  N  u 

ill  Xt  +  b  i£i  Xt  +  c  !«!  XjYj  -  i=!  X1Z1 


«  E  +  b  z  X.Y<  +  c  j  yJ  •  J,  Y.Z4 
J~'  1  i»l  1  1  i»l  1  1*1  1  1 


N 

iSi 


(17) 


H  N  N 

aN  +  b  1E1  Xi+  c  ^  Yt  -  Z± 


Local  coordinate  systea 

61.  Now  that  a,  b,  and  c  are  known,  the  new  set  of  points  for  the 
modified  base  (plane  of  analysis)  is  computed  froa  Equation  16.  It  should  be 
noted  that,  if  the  original  points  lie  in  a  plane,  the  new  points  will  be  the 
aaae  as  the  old  ones.  Next,  the  global  coordinates,  forces,  and  moments  are 
transformed  to  a  local  coordinate  systea  where  X^  and  Y^  lie  in  the  plane  of 
analysis  and  Z^  is  normal  to  the  plane  of  analysis.  This  transformation  is 
made  by  two  rotations  and  a  translation.  First,  as  shown  in  Figures  9  and  17, 
a  rotation  is  done  about  the  Z  axis  so  that  the  new  X'  axis  is  parallel  to  the 

downwardmost  direction  of  the  modified  base. 


V 


Figure  17.  Z  axis  rotation 


Mathematically, 

[X*  Y'  Z’]  ■  (X  Y  Z] 


cosd  „ 

“Sin0.„ 

o' 

xy 

xy 

sinO 

cos©_ 

0 

*y 

*y 

0 

0 

1 

(18) 


i 

Next,  as  shown  in  Figures  9  and  18,  a  Y*  axis  rotation  is  done,  as  follows: 


N 


IX"  Y"  Z"]  -  [X*  Y'  Z’] 


CO80 


0  sin© 


zn  ~  zn 

0  10 

-sin0„  0  cos9 

zn  zn 


Figure  18.  Y'  axis  rotation 
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Last,  as  shown  in  Figures  9  and  19,  a  translation  to  have  the  plane  of 

ZL 


Figure  19.  Base  in  local  system 
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62.  Uplift  pressures  can  be  specified  at  the  points  defining  the  base  of 
investigation.  When  the  original  points  are  replaced  by  a  different  set  of 
points  that  lie  in  the  plane  of  analysis,  the  uplift  values  must  also  be 
modified.  The  equation  used  is 

Unew  "  Uold  +  Yw  <Zold  “  Znew>  <21> 

where  is  the  density  of  water. 

63.  Next,  the  contribution  of  uplift  to  the  other  forces  and  moments 
must  be  computed.  In  addition  to  the  uplift  pressures  being  specified  at 
certain  points  on  the  boundary,  some  assumptions  must  be  made  as  to  what 
values  the  uplift  has  at  other  points  on  the  base. 

64.  First,  uplift  is  assumed  to  vary  linearly  along  the  line  segments 
between  the  specified  points.  Next,  the  base  is  rotated  (for  computation 
purposes  only)  so  that  the  line  segment  with  the  maximum  gradient 


L 


is  parallel  to  the  local  X  axis.  Here,  and  Uj  are  the  uplift  values 
specified  at  the  two  ends  of  a  line  segment,  and  L  is  the  length  of  the  line 
segment.  Next,  as  shown  in  Figure  20,  the  base  is  divided  into  horizontal 
strips.  Note  that  a  horizontal  line  is  always  created  from  a  vertex  and 


V 


sometimes  midway  between  vertices.  Note  also  that  vertical  lines  are  drawn 
when  the  midlines  are  needed  to  further  subdivide  the  base.  Uplift  pressures 
are  assumed  to  vary  linearly  along  these  horizontal  lines. 


v 

Line  Segment  with  Maximum  Sradient 

Figure  20.  General  base 

Combining  that  assumption  with  :he  linear  variation  along  the  boundary 
assumption  allows  all  the  uplift  values  to  be  computed  at  the  vertices  of  each 
of  the  constructed  triangles  or  quadrilaterals.  For  Instance,  is  computed 
from  U^  and  Uj,,  and  is  computed  from  Ug  and  Uj.  Then,  can  be  computed 
from  Ug  and  U^. 
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TVt’v'  1MJI 


compute  the  uplift  for  any  point  inside  a  triangular  element,  a 
:ion  is  assumed.  That  is. 


U  -  a  +  bXL  +  cYl 

i»nd  c  are  constants  that  can  be  evaluated  from  the  corner  values, 
r  the  quadrilateral  element,  the  isoparametric  variation  is 

U  »  |  <1  -  S)(l  -  T)U1  +  |(1  +  S) (1  -  T)U2 

+  \  (1  +  S)(l  +  T)U3  +  |(1  -  S)  (1  +  T)U4 


,  U3,  and  U4  are  the  corner  values,  and  S  and  T  are  parameters 
stween  -1  and  1.  For  the  case  of  a  rectangular  base  where  =  U4 
,  the  above  equation  reduces  to: 


U  =  -|  (1  -  S)UX  +  \  (1  +  S)  u2 


+  u2)  + 


S 


mple  linear  interpolation. 

inally,  to  compute  the  contribution  of  uplift  to  the  forces  and 
is  merely  necessary  to  integrate  over  each  of  the  individual 
Since  uplift  acts  normal  to  the  base,  only  a  Z  component  of  forces 
he  local  coordinate  system.  Thus, 

Fu  “  /AUda  (22) 

ent  contributions  are 

Mxu  “  4YLUda 

“yu  "  -Wda  <23> 
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APPENDIX  A:  BUILDING  AND  MODIFYING  DATA  FILES 


1.  This  appendix  contains  some  general  instructions  on  how  to  build  and 
modify  data  files. 

WES  and  Macon  Computer  Systems 

2.  To  build  a  new  data  file,  the  user  first  types 

NEW 

He  then  types  the  file  with  line  numbers  until  finished.  The  *  will  appear 
after  each  line  is  typed.  Next,  once  the  file  is  completed, 

SAVE  NAME 

is  typed,  where  NAME  is  a  file  description.  If  a  line  is  wrong,  the  user 
simply  retypes  it.  Then,  he  types 

RE SAVE  NAME 

after  making  all  corrections. 

3.  Input  data  files  for  this  program  are  read  with  a  free-field 
format.  Thus,  data  are  separated  by  either  blanks  or  commas,  and  floating¬ 
point  data  that  are  whole  numbers  do  not  need  a  trailing  decimal  point.  The 
following  is  a  sample  data  file  building  sequence  In  the  style  of  the  U.  S. 
Arm:  Engineer  Waterways  Experiment  Station  (WES)  and  the  Office  of  Personnel 
Management,  Macon,  Ga.  (Macon),  computer  systems: 


NEU 

X10  1  10.  0. 

*20  2  10.  10. 
*30  30.  10 
*40  0.  0. 

XSAUE  DAM1 
DATA  SAUED-DAM 
*5  POINTS  4 
XRESA-  DAM 
DATA  SAUED-DAM 
XLIST 


S  POINTS  4 
10  1  10.  0. 
20  2  10  10. 
30  3  0.  10. 
40  0.  0. 


Control  Data  Corporation  (CPC)  Computer  System 


4.  To  build  a  new  data  file,  the  user  first  types 


NEW, NAME 


where  NAME  Is  a  one-  to  seven-character  name  starting  with  a  letter.  Then,  he 


types  the  file  with  line  numbers  until  finished.  The  C>  prompt  will  not 


reappear  after  the  first  line  is  typed.  The  user  next  types 


SAVE, NAME 


upon  completion.  If  a  line  is  wrong,  he  simply  retypes  It.  Finally,  he  types 


REPLACE, NAME 


after  making  all  corrections. 


5.  Input  data  files  for  this  program  are  read  with  a  free-field 


format.  Thus,  data  are  separated  by  either  blanks  or  commas,  and  floating¬ 


point  data  that  are  whole  numbers  do  not  need  a  trailing  decimal  point.  A  CDC 


computer  system  sample  data  file  building  sequence  appears  as  follows: 


A2 


C>NEU,L0CK1 
C>10  POINTS  4 
20  1  0  0 
30  2  100.5  0 
40  3  100.5  20 
50 

SAUE>L0CK1 
OLIST 
10  POINTS  4 
20  1  0  0 
30  2  100.5  0 
40  3  100.5  20 
050  4  0  20 
REPLACE, LOCK 1 
OLIST 
10  POINTS  4 
20  1  0  0 
30  2  100.5  0 
40  3  100.5  20 
50  4  0  20 
C> 


A3 


In  accordance  with  letter  fro*  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsinile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Tracy,  Fred  T. 

A  three-dimensional  stability  analysis/design  program 
(3DSAD)  :  Report  3  :  General  Analysis  Module  (CGAM)  /  by 
Fred  T.  Tracy  (Automatic  Data  Processing  Center,  U.S. 

Army  Engineer  Waterways  Experiment  Station)  and  Charles 
W.  Kling  (U.S.  Army  Engineer  District,  Mobile).  —  Vicksburg, 
Miss.  :  The  Station  ;  Springfield,  Va.  :  available 
from  NTIS,  1982. 

51,  3  P.  :  ill.  ;  27  cm.  —  (Instruction  report  ;  K-80-4, 
Report  3) 

Cover  title. 

"June  1982." 

"A  report  under  the  Computer-Aided  Structural 
Engineering  (CASE)  Project." 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army." 

Bibliography:  p.  SI. 

1.  Computer  programs.  2.  Stability.  3.  Structural 
stability.  4.  3DSAD  (Computer  program).  I.  Kling, 

Charles  W.  II.  United  States.  Army.  Corps  of  Engineers. 
Office  of  the  Chief  of  Engineers.  III.  U.S.  Army 


Tracy,  Fred  T. 

A  three-dimensional  stability  analysis/design  :  ...  1982. 

(Card  2) 

Engineer  Waterways  Experiment  Station.  Automatic  Data 
Processing  Center.  IV.  United  States.  Army.  Corps  of 
Engineers.  Mobile  District.  V.  Title  VI.  Series: 
Instruction  report  (U.S.  Army  Engineer  Waterways  Experiment 
Station)  ;  K-80-4,  Report  3. 
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WATERWAYS  EXPERIMENT  STATION  REPORTS 
PUBLISHED  UNDER  THE  COMPUTER-AIDED 
STRUCTURAL  ENGINEERING  (CASE)  PROJECT 


Technical  Report  K-78-1 
Instruction  Report  0-79-2 

Technical  Report  K-80-1 
Technical  Report  K-80-2 

Instruction  Report  K-80-1 

Instruction  Report  K-80-3 
Instruction  Report  K-80-4 

Instruction  Report  K-80-6 
Instruction  Report  K-80-7 
Technical  Report  K-80-4 

Technical  Report  K-80-5 
Instruction  Report  K-81-2 

Instruction  Report  K-81-3 
Instruction  Report  K-81-4 
Instruction  Report  K-81-6 

Instruction  Report  K-81-7 
Instruction  Report  K-81-9 
Technical  Report  K-81-2 
Instruction  Report  K-82-8 


Title 

List  of  Computer  Programs  for  Computer-Aided  Structural  Engineering 

User’s  Guide:  Computer  Program  with  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

Survey  of  Bridge-Oriented  Design  Software 

Evaluation  of  Computer  Programs  for  the  Design/Analysis  of 
Highway  and  Railway  Bridges 

User’s  Guide:  Computer  Program  for  Design/Review  of  Curvi¬ 
linear  Conduits/Culverts  (CURCON) 

A  Three-Dimensional  Finite  Element  Data  Edit  Program 

A  Three-Dimensional  Stability  Analysis/Design  Program  (3DSAD) 
Report  1:  General  Geometry  Module 
Report  3:  General  Analysis  Module  (CGAM) 

Basic  User's  Guide:  Computer  Program  for  Design  and  Analysis 
of  Inverted-T  Retaining  Wails  and  Floodwalls  (TWDA) 

User's  Reference  Manual:  Computer  Program  for  Design  and 
Analysis  of  Inverted-T  Retaining  Walls  and  Floodwalls  (TWDA) 

Documentation  of  Finite  Element  Analyses 
Report  1:  Longview  Outlet  Works  Conduit 
Report  2:  Anchored  Wall  Monolith,  Bay  Springs  Lock 

Basic  Pile  Group  Behavior 

User’s  Guide:  Computer  Program  for  Design  and  Analysis  of  Sheet 
Pile  Walls  by  Classical  Methods(CSHTWAL) 

Report  1:  Computational  Processes 
Report  2:  Interactive  Graphics  Options 

Validation  Report:  Computer  Program  for  Design  and  Analysis  of 
Inverted-T  Retaining  Walls  and  Floodwalls  (TWDA) 

User's  Guide:  Computer  Program  for  Design  and  Analysis  of 
Cast-in-Place  Tunnel  Linings  (NEWTUN) 

User’s  Guide:  Computer  Program  for  Optimum  Nonlinear  Dynamic 
Design  of  Reinforced  Concrete  Slabs  Under  Blast  Loading 
(CBARCS) 

User’s  Guide:  Computer  Program  for  Design  or  Investigation  of 
Orthogonal  Culverts  (CORTCUL) 

User’s  Guide:  Computer  Program  for  Three-Dimensional  Analysis 
of  Building  Systems  (CTABS80) 

Theoretical  Basis  for  CTABS80:  A  Computer  Program  for 
Three-Dimensional  Analysis  of  Building  Systems 

User's  Guide:  Computer  Program  for  Analysis  of  Beam-Column 
Structures  with  Nonlinear  Supports  (CBEAMC) 


Date 
Feb  1978 
Mar  1979 

Jan  1980 
Jan  1980 

Feb  1980 

Mar  1980 

Jun  1980 
Jun  1982 

Dec  1980 
Dec  I960 

Dec  1900 
Dec  1980 

Dec  1980 

Feb  1961 
Mar  1961 

Feb  1961 
Mar  1961 
Mar  1981 

Mar  1961 
Aug  1961 
Sep  1961 
Jun  1962 


